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Abstract. The following paper is a variation on a theme of Gian- 
franco Cimmino on some integral representation formulas for the 
solution of a linear equations system. 

Cimmino was probably motivated for giving a representation 
formula suitable not only for theoretical investigations but also for 
applied computation. 

In this paper we will prove that the Cimmino integrals are 
strictly related to the residues of some zeta-like functions asso- 
ciated to the linear system. 



1. Introduction 

Gianfranco Cimmino was born in Naples on March 12, 1908. 

He received his Laurea degree in Mathematics at the University of 
Naples under the direction of Mauro Picone (1885-1977) in 1927. 

At the end of 1939 Cimmino moved permanently to the University 
of Bologna to occupy the chair of Mathematical Analysis. 

Cimmino died in Bologna on May 30, 1989. 

Gianfranco Cimmino (1908-1989) was a student of Renato Cacciop- 
poli (1904-1959) and Mauro Picone toghether with Giuseppe Scorza 
Dragoni (1907-1996) and Carlo Miranda (1912-1982). 

Mauro Picone founded the "Istituto per le Applicazioni del Calcolo" 
(IAC) in 1927. Indeed, long before the introduction of digital com- 
puters, Picone had the intuition of the potential impact on real-life 
problems of the combination of computational methods with mathe- 
matical abstraction. 

Probably the influence of Picone ideas and Cesari papers [2J and [3J 
on numerical solution of linear systems leads Cimmino to be interested 
to the numerical treatment of the solutions of linear systems of algebraic 
equations. 

Among his researches he also obtained an interesting representation 
of the solution of a linear system of real linear equations that now we 
describe. 

Let 



(1) 



Ax = b 
l 



be a linear system of n equation and n unknown, where the unknown 
values the components if the column vector x G JR n and 

A is non singular matrix with real coefficient of order n. 

The well known Cramer's rule say that if D = det(A) ^ then 

(2) Xi = ^j, i = l,...n, 

where D{ is the determinant of the matrix obtained replacing the %— th 
column of the matrix A with the column vector b. 

The alternative representation of ([2]) given by Cimmino is 

(3) Xi = i = l,...n, 
where 

(4) C = [ \\A f uf n du, 
and 

(5) Ci = n \\A t u\\~ n ~ 2 (b, u) (A f x, e*) dix. 

The integration here is made with respect to the standard (n — 
1)— dimensional measure on S 1 ™ -1 , the boundary of the Euclidean unit 
ball in El", ||-|| and (•, •) stand respectively for the Euclidean norm and 
the Euclidean inner product on H n and e\, . . . , e n is the canonical basis 
of H n , and A 1 denote the transpose of the matrix A. 

In [4j, [5] and |7] Cimmino gives a probabilistic argument which 
justify the existence of such kind of formulas and in |S] he also give 
an elementary but not trivial proof of See also ([S]) for some 

applications. We refer to [1] for a detailed discussion and a background 
information on Cimmino's papers in the field of the numerical analysis. 

The purpose of this paper is to show that the Cimmino ideas fit 
nicely into the theory of the residues of zeta-like functions, a powerful 
tool coming from the analytic number theory. 

Namely we will show that (jl]) and (jSj) are the integral representation 
of the residues of suitable zeta-like function associate to the matrix 
A and the vector b of the linear system ([T]); see Theorem 17.21 for the 
complete statement. 

It should be noted that formulas H] and [5] actually are due to Jacobi: 
cf. [TO]. 

2. Theta functions and their Mellin transform 

Let us begin with the following (almost trivial) observation: 

Proposition 2.1. Let f be a continuous complex function defined on 
the real interval [0,1]. Assume that for some constants R,a,(3 G 1R 
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with a > f3 we have 

fit) = Rt~ a + 0(t~ p ), t^0 + . 
Then, given sG C, the integral 

9(s) = jf /(*)*• j 

converges absolutely on the half space Re(s) > a and extends to a 
meromorphic function on the half space Re(s) > (3 having a simple 
pole at s = a with residue R. 

Proof. Since f(t) = 0(t~ a ) then g(s) converges and is holomorphic 
when Re(s) > a. If we denote 

fc(t) = f{t)-Rr a 

then 

,( S ) = / ,1 /(^y = — + /Wy- 
Jo t s-a Jo t 

Since = 0(t~' 3 ) the last integral define a holomorphic function on 
the half space Re(s) > /3 and we are done. 
□ 

Thus the singularities of the analytic function g(s) describe the be- 
haviour of the function fit) as t — > + . 

Let us recall that C(R") and L 1 (IR n ) denotes respectively the space of 
the continuous complex function on 1R™ and the space of the absolutely 
integrable complex functions with respect to the Lebesgue measure on 
R n 

Given / e L 1 (R n ) the Fourier transform of / is the function / 
defined by the formula 

f(y) = [ f(x)e- 2 ^>dx. 

Observe that 

/(0) = / f(x)dx. 

We denote by 5(lR n ) is the Schwartz space of smooth functions 
rapidly decreasing at the infinity together with their derivatives of all 
orders. Given a, (3 > two positive real number we denote by 

s:t" a (R n ) 

the space of all measurable function such that 

sup \f(x)\(l + \\x\\ n+a ) <+oo 
xeR" 

and 

sup 

3 



f(y) (1 + \\y\ n+p ) < + 



oo. 



Of course (R n ) C L l (R n ) and S(R n ) is the intersection of all the 
spaces «S^f(E n ) when a and /3 varies on all the positive real numbers. 
Since 

a„3>0 



is useful to set 



Let a, /3 > 0, finite or infinite, be fixed . 
Let / G (M n ) be an arbitrary function. 

By the very elementary approach to the (Riemann) integration the- 
ory we have 

/(0)= / f(x)dx = lim f Yl /M- 
It is then natural to define for t > and d > 

W,*) = E 'f(t 1/d ") = 8<i (/,*)- /(0), 
wez n 

(where XLez™' as usua l ; stands for XLez™\{o}) an d considering the 
Mellin transform. 

f +oc dt 

= J e*Af,t)t° T 

(The introduction of the constant d will simplify some computations). 

Our hope is to study the Riemann approximation looking at the 
residues of £ d (/, s). 

Indeed we have 

Theorem 2.1. Let a, f3 > be given and let f E S£%(R n ) be an 
arbitrary function. Then the integral defining ^ (/, s) converges in the 
strip n/d < Re(s) < (n + a)/d and the function £^ (/; s) extends to a 
meromorphic function on the strip —(3/d < Re(s) < (n + a)/d. having 
exactly two simple poles respectively at s = n/d with residue /(0) and 
at s = with residue — /(0). Moreover we have the functional equation 

(3 . . n + a 



& (f, s) = £ d (J, n/d - s) , --< Re(s) < 



d 



The proof of the theorem above follows closely the lines of (some 
of) the standard proof of the functional equation for the classical zeta 
functions used in number theory, but the statement of the theorem in 
the form above is not so common. We refer [IT] for a detailed proof. 
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Observe that 6* d (/, t) = 6\ (/, t 1 ^) and hence, by a simple change of 
integration variable, ^ (/, s) = d(,i (/, ds). It follows that the function 
£d (/, s) has a pole at s = n/d if, and only if, the function £1 (/, s) has 
a pole at s = n with the same residue. 

3. Zeta functions machinery 

The following proposition describe the hearth of our approach to the 
treatment of Cimmino integrals. 

Proposition 3.1. Let a,(3,d > be given and let f E S(R n ) be an 
arbitrary function. Assume that in the strip —j3/d < Re(s) < (n+a)/d 
the function ^ (/, s) decomposes as 

= Gf(s)Z f (s) 

with Zf(s) being meromorphic with exactly a single simple pole at s = 
n/d. 
Then 

Z f (0)ResG f (s) = -/(0), 

s=0 

G f (n/d) Res Z f (s) = /(0). 

s=n/a 

Proof. By the previous theorem £ d (/, s) has a simple pole at s = and 
-/(0) = Res6(/, S ) = ResG / ( S )^/( S ). 

Since Zf(s) is holomorphic at s = then 

-/(0) = Z f (0)ResG f (s). 

s=0 

Replacing / with /, we have 

-/(0) = Res6(/,s). 

Using the functional equation t; d (/, s) = £ d n/d — s^j we obtain 

/(0) = - Res ^ (f, s)=- Res ^ (/, n/d - s) = Res & (/, s) . 

Since ^ (/, s) and -Z/(s) have a simple pole at s = n/d then necessarily 
Gf(s) is holomorphic at s = n/d and hence 

/(0) = G/(n/rf) Rffl//(«), 

as desired. 
□ 

In the sequel we will apply the above proposition to functions / which 
gives a decompositions ^ (/, s) = Gf(s)Zf(s) where the function Gf(s) 
is an algebraic combination of elementary functions and functions of 

5 



the form T(as + b), being T(z) the Euler gamma function, and Zf(s) 
is representable when Re(s) >> as (generalized) Dirichlet series 

oo 

Zf(s) = $>A- 

k=l 

where Cfc, G IR, Ai < A2 < • • • , and A& — > +00 as k — > +00. 

4. Quadratic forms 

We begin considering the well known classical example of the Gauss- 
ian integrals associated to quadratic form. 

Let us recall that T(s) denotes the Euler gamma function; it is a 
meromorphic function on C which for Re(s) > satisfies 

r + °° dt 

!■(.) = / A-* 

We assume the reader knows all the properties of such a function. 
Let Q be a real symmetric matrix of order n. We set 

qq{x) = (Qx,x) 
Assume that Q is positive definite and consider the function 

g Q (x) = e-^W 
then gq E S(M. n ) and a standard argument yields 

In particular 

Since qg(tx) = t 2 qq(x) then, choosing d = 2, 

wez n 

and hence, when Re(s) > n/2, 

wez n ^° 

The change of variable u = tnqqiuj) yields 

/ e ~ UuS - E^h _s 



6 



The function defined for Re(s) > n/2 by 

is the Epstein zeta function associated to the quadratic form qg. 

Since gg{x) G iS(IR n ) then, by Theorem 12. 1[ £2 {9q,s) is a meromor- 
phic function on (D having exactly two simple poles at s = n/2 and 
s = with residues respectively R = 1 and R = —1. 

Since the function 7r _ T(s) never vanishes on (D and the function T(s) 
has a simple pole at s = with residue R = 1 it follows that £ (gg, s) 
extends to a meromorphic function on (D having exactly a simple pole 
at s — n/2. 

By Proposition 13.11 and the relation 

r(s + l) = sT(s) 
we obtain that the residue of ( (qg, s) at s = n/2 is 

vr"/ 2 1 n ( 7r n / 2 \ 1 



Res t(q Q ,s) 



s=n/2 "w > r(f)Vdetg 2 \^r(f + 1)/ v/deTQ 
and that 

C(QQ,s) = -1. 
All that is well known to number theorists. 

But now consider two quadratic forms, qg and q$ with Q and 5* 
symmetric matrices with Q positive definite and set 

9qAx) = qs(x)e~^ x) . 

Then g Q;S E S(R n ). 

Proposition 4.1. The Fourier transform of gg y s{x) is 
* , s 1 r ^ Tr(<5 _1 ^) / x 

1 m-^Q^e-^-^ + ^^^e'^\ 



VdetQ v w/ 2vr v ^etQ 

where Tt(Q~ 1 S) denotes the trace of the matrix Q~ 1 S . 
In particular we have 

Tr(Q~ 1 S) 

gg,s{x) dx = g Q , s (0) = - — j^=^- 



— (x, x) — ^ ] x i , 



Proof. The function 
where 

II II 2 

If II 

i=i 

coincide with its Fourier transform. Set 

Then, the standard properties of the Fourier transform yields 

where D y . denote the operator of derivation with respect to the variable 
Vi- 

We compute 

My) = ViV^ AvX \ 
D v Mv) = -^vfi Avl \ 

D yi D yj h tJ (y) = -27rD yi ( yj e-^ 2 ) = 47t 2 e~^ 2 - 2n5 tJ e^\ 
My) = —Dy.D^M 2 = -e^ + 2_ Sije ^\\ 

Let denote by I n the identity matrix of order n. Let C = (c^) be a 
real symmetric matrix of order n and set 

n 

hc(x) = gi n ,c( x ) = qc(x)e-* M2 = ^ c l3 h l3 [x) 



Then 



hc(y) = ^Cijhijiy) 



^ ] C ij e ^ ] r, C ij^ij e 



2tt 

ij=i ij=i 
-gc(y)e-^ l|2 + ^Tr(C)e- 



Let Qz be the unique positive definite symmetric matrix sue that 

{q^) 2 = q 



and let Q 2 be its inverse. 
Set 



f(x) = h c (x) 



1 ^ ^ 1 



where C — Q 2 SQ 2 . We have 

Tr(C) = Tr(g-^g-i) = Tr(Q-§g-^) = Tr^S) 

and hence 



/(v) = -qc{y)e^ 2 + ^{Q- 1 s)e-^ 2 . 

In 



We also have 



Qc{Q 2 x) = qs(x), 

qc(Q~*y) = ?Q-i5Q-i(y), 



?<?(2/)- 



It follows that 

and hence 
9Q,s(v) 



9q,s( x ) = f(Q*x) 



-J(Q-*y) 



v/detg- 

1 -^s Q Ay)e-^ {y) + ^Me~^ 



v/detg 



27r^/detQ 



as desired. 
□ 



Let now compute 9* (g Q>s , t) and £ 2 (g Q>s , s). 

Since qqitx) = t 2 qq(x) and qsitx) = t 2 qs(x) we have 



e *2 (dQ,s,t) = t( ls(uj)e- 



tirq s (ui) 



and 



r+co 1 

J° ,,,C7« 

= E 



dt 
T 



+00 7, 
qsi^e-^^h^ 1 -. 
* 



By the change of variable u = tnqs(uj) we obtain 



(r+oo 



e~ u u s+1 



dt 



qs(oo)q Q (co)-^ 



Set 



/ 

((<1q,<1s,s) = Y Qs(u)qq(w)' 



wez n 
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The meromorphic function 7r _ ^ +1 ^r(s+l) never vanishes on (D. Since 
9q,s{Q) — by Theorem 12.11 the function £2 {9q,s,s) is a meromorphic 
function on (D having exactly a simple pole at s = n/2 and hence the 
sum defining ( (qQ,qs, s) converges for Re(s) > n/2 and extends to a 
meromorphic function on (D having exactly a simple pole at s = n/2. 
Then proposition 13. II and proposition 14.11 give: 

Theorem 4.1. Let Q and S be two symmetric matrices of order n with 
Q positive definite. Then the sum defining the zeta function ( (qQ, q$, s) 
converges absolutely for Re(s) > n/2 + 1. 

The zeta function ( (qQ, q$, s) extends to a meromorphic function on 
(D having exactly a simple pole at s = n/2 + 1 with residue 



Tr{Q- l S) 




J$H C ^'™'~' 2\T(f + l)y v/detQ 
From the equality 

9qM = -7m gQ ^ Q ~ 1SQ ~ l{y) + 2^jm gqAy) 

proved in Proposition 14. II we obtain 

&(9q,s,s) = --^=n-^r(s + l)c{q qQ ^,q Q -r S Q-^s + l 



Inserting such expressions in the functional equation 

6 (9q,s, 7; - s ) =6 (9q,s, s) 
and dividing by it we obtain: 

Theorem 4.2. Let Q and S be as in theorem \4-l\ The zeta function 
C {qQ) qsi s ) satisfies the functional equation 

'n \ , / n 



n -K2-')T + 1 " *J C {QQ> qs,^ + l-s 

ir- s r(s + i)((q Q -i,qQ-i S Q-i,s + l) 



VdetQ 
^ Tr(Q-^) 

2 x /detQ V A Wg ; 

5. Lattices 

A lattice in IR n is a set of the form 

L = {Auo I w 6 Z n } 

where A e GL(n, 1R) is a real invertible matrix of order n. 

If A x is an other matrix such that Ai(Z n ) = L then U = AA^ l (Z n ) C 
Z n . It follows that the invertible matrix U has integral coefficients and 
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hence |det U\ = 1, that is |det A| = |detAi|. We define the volume of 
the lattice L as 

\L\ = \detA\ . 

By the argument given above the definition does not depend on the 
choice of the matrix A. 

If L = A(Z n ) is a lattice the dual lattice is the lattice L A associate 
to the inverse of the transpose of the matrix A. For convenience we 
also set 

A = (A 1 )' 1 

Given a lattice L C !R n and a positive definite symmetric matrix Q 
of order n we define 

(L{q Q ,s) = 9qM~* 
and if S is any symmetric matrix we also define 
Cl (gg, qs, s) = <ls(u)qQ(uy s - 
Of course we have 

Cz" (qQ,s) = ((q Q ,s) 

and 

Cz™ (qQ,qs,s) = C(QQ,Qs,s) ■ 
If L = A(Z n ) with A e GL(n, JR) then we have 

Cl (Qq, s) = C (q^QA, s) , 
where A 1 denotes the transpose of the matrix A, and 

Cl (qq, qs, s) = C {qA*QA, qA*sA, s) . 

We also have 

y/de^QA) = |det A\ ^detQ = \L\ VdetA 

and 

Ti({A t QA)-\A t SA){A t QA)- 1 ) = Tr (A' 1 (Q-'SQ' 1 ) A) 

= TviQ^SQ- 1 ). 
Applying the theta-zeta machinery to the function 

9a,q,s(x) = 9Q,s{Ax) = qs(Ax)e-^ Ax \ 
observing that the Fourier transform of gA,Q,s( x ) is 

9a,q,s(v) = | de | A] 9Q,s{Ay), 

it follows that the results of the previous section generalize: 
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Theorem 5.1. Let L C IR™ be a lattice and Let Q and S be two sym- 
metric matrices of order n with Q positive definite. Then the sum 
defining the zeta functions (qq, s) and (l {qQ, q$, s) converges abso- 
lutely respectively for Re(s) > n/2 and Re(s) > n/2 + 1. 

The zeta functions Ql {qQ, s) and Ql {qQ, qs, s) extends to a meromor- 
phic function on (D having exactly a simple pole respectively at s = n/2 
and s = n/2 + 1 with residues 



s n I Ti 

Res Cl {qQ, s 



n/2 



and 



•=n/2^ v ^' ' 2 \ T{± + 1)J \L\ v / detQ' 

. , \ Trrg- 1 ^) 

Res Cl {qQ, qs, s) = - . -r /xTTr 
^=«/2+i 2 I T (f + 1) / \L\ VdetQ 



Theorem 5.2. Lei L C H n , Q and S be as in theorem \5.l\ The zeta 
functions Cl{Qq,s) and Cl (O'Q; <?s> s ) satisfy the functional equations 
respectively 



71 



1 

:7T T(s)CiA (g Q -l,s) 



and 



\L\ VdetQ 



7r -(f-)r(| + l-s)a(gg,gs,| + l 



7T T(s + 1)Cla (g Q -i,g Q -i SQ -i,s + l) 



\L\ y/detQ 
Tr{Q^S) 



2\L\ v^ietQ 



ti T(s)Cla (9q-i,s) . 



6. Integral representation 

Let gQ and qs be two quadratic form with Q and S 1 symmetric ma- 
trices and Q positive definite. 

We now will give such residues as integrals over the boundary of the 
unit ball in IT. 

We already observed that 

e-^ x) dx 1 



/ R n y/detQ 
and by Proposition 14.11 we also have 

R n 27rvdet Q 
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If A e GL(n, IR) then a simple change of variable gives 

1 



e 



-nq Q (Ax) j x 



V |det A\ y/det Q 

and 



Qs 



(Ax)e~^ {Ax) dx 



/ R » 2n |det A\ y/det Q 

Let us recall that given / 6 L 1 (lR n ) the integration by polar coordi- 
nates gives 

f(x)dx= / I / f (ru)^" 1 dr ) du, 
where dtt is the Euclidean (hyper-) surface measure on the unit sphere 

s n -\ 

Using such formula we obtain 

Tt (Q~ 1S ) = f „JAx)e-^ Ax) dx 

27r|detA| y/dKQ Jn- 



J Qf q s (rAu)e-^ {rAu) r n - l dr^j du 



qs(Au) I J e -^u) r n + 2±}j du 

By the substitution nr 2 qQ(Au) = t in the inner integral we obtain 

[ q s {Au)( ! + °° e-^ Au) r n+2< ^-\ du 
Jsn-i \Jo r J 



and hence 



g 5 (^)gQ(^)" (n/2+1) du 

gn-l 



71 



n/2 



r(f + 1) J \detA\ y/deTQ' 



When S = Q we obtain 

q Q (Au)~ n/2 du = n 



r (2 + 1) y |det A| v/Hetg' 

Replacing s with s/2 in the formulas of the residues in Theorem 15.11 
we obtain 



ResCL(<?Q,s/2) = n 



7T 



n/2 



r(f + 1)1 \L\ v / deTQ ! 
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and 



Res Cl (qQ,qs,s/2) 



71 



n/2 



S =n+2^ V ^'^' ' > \Y{\ + l)J \L\ v^tQ' 

Thus we obtained: 

Theorem 6.1. Let qq and q$ be two quadratic form with Q and S 
symmetric matrices and Q positive definite. Let A 6 GL(n, IR) and set 
L = A(Z n ). Then we have 

Resa(<Zg,s/2)= / q Q (Au)- n / 2 du. 

s=n J gn-l 



and 



Res C L (q Q ,qs,s/2)= [ q s (Au)q Q (Au)- (n / 2+1 Uu, 

s=n+2 JS 71 - 1 



7. Linear systems 
We recall that A denotes the inverse of the transpose of the matrix 

A. 

Given A e GL(JR n ) we define 



((A,s)= \\M\~ 



U)£Z n 

and given also b G IR n we define the vector value zeta function 
C(A,b,s)= J2'\\M\~ 2s (b^)Au. 

Of course we have 

((A, s) = ( L (q In ,s) 

with L = A(Z n ), I n the identity matrix of order n and hence, by Theo- 
rem [5J3 we have 

t- 1 a \ n ( 7t n/2 \ 1 

Res C(A,s) = -\ —7 r tz j7 

s=n/2 SK ' 2 \ T(f + l)J \detA\ 

But we also have: 

Theorem 7.1. The series defining ( (A, b, s) converges for Re(s) > 
n/2 + 1 and the function ((A,b,s) extends to a (vector value) mero- 
morphic function on (D having only a simple pole at s = n/2 + 1 with 
residue 

1 / W 2 \ 1 * 

Res C (A, b,s) = -\ —. 7 -Ab. 

s=n/2+r K ' ; 2 \r(f + 1)J \detA\ 
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Moreover, if c G IR n then we have the functional equation 
n-G-frfe + l-s) (c(A,b^ + l-s),Ac 



2|det A\ v 7 V ' J 

Proof. Given u, v G IR n , with w = (%,... , tt„) and v = (v±, . . . , v n ) 
we denote by u <8> f the symmetric matrix (?%•) of order n with entries 
given by Wij = UiVj. Then we have 

Tr(w g) u) = (u, v) . 

Now observe that 

(((A,b,s),Ac) = ( L (q In ,q s ,s) 
where L = A(Z n ), I n is the identity matrix of order n and 

S = Ab ® Ac. 

We have 

Tr(ife <g> Ac) = Ac\ = (b, c) 



and hence, the formulas for the residue and the functional equation 
of the function (A, 6, | + 1 — s) , Ac) follows from Theorem 15.11 and 
Theorem 15.21 
□ 

We are now ready to state and prove the reformulation of Cimmino's 
results. 

Theorem 7.2. Let 

Ax = b 

be a linear system of n equation and n unknown, where the unknown 
values the components if the column vector x G IR n and 

A is non singular matrix with real coefficient of order n. 
Then we have 

R 

X{ 

where 



Xi , % 1, ... 71, 



R = ResC(A*,s/2) 

s=n 

and 

Ri = n Res <( (A\ b, s/2) , a) , 

s=n+2 

where e\, . . . , e n is the canonical basis o/]R n . 

15 



Moreover we have the identities 

R = [ HA*ulP' du, 

and 

Ri = n ||A*w||~ n ~ 2 (b, u) (A t x,e i ) du. 
Proof. By assumption 

x = A~ x b 

and we have 

where ei, . . . , e„ is the canonical basis of IR™. As previously observed 
we have 

1 




|det A\ 



and hence 



R 

as required. 

The last assertions of the Theorem follows from Theorem 16.11 
□ 

Our tour around Cimmino's ideas is completed. 
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